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Particles in space periodic potentials constitute standard models for investigation of crystalline 
phenomena in solid state physics. Time periodicity of periodically driven systems is a close analogue 
of space periodicity of solid state crystals. There is an intriguing question if solid state phenomena 
can be observed in the time domain. Here we show that wave-packets localized on resonant classical 
trajectories of periodically driven systems are ideal elements to realize Anderson localization or Mott 
insulator phase in the time domain. Uniform superpositions of the wave-packets form stationary 
states of a periodically driven particle. However, an additional perturbation that fluctuates in time 
results in disorder in time and Anderson localization effects emerge. Switching to many-particle 
systems we observe that depending on how strong particle interactions are, stationary states can 
be Bose-Einstein condensates or single Fock states where dehnite numbers of particles occupy the 
periodically evolving wave-packets. Our study shows that non-trivial crystal-like phenomena can be 
observed in the time domain. 


Nearly sixty years ago Anderson discovered that trans¬ 
port of non-interacting particles in the presence of disor¬ 
der can stop totally due to interference of different scat¬ 
tering paths [l|. The suppression of the transport is ac¬ 
companied by localization of eigenstates in the configu¬ 
ration space ii- Quantum description of classically 
chaotic systems reveals yet another version of Anderson 
localization (AL) where localization of a particle takes 
place in the momentum space [3|. Suppression of classi¬ 
cal chaotic diffusion in the momentum space was identi¬ 
fied with AL after the quantum kicked rotor system was 
mapped to the quasi-random one-dimensional (ID) An¬ 
derson model . Here we show that the presence of time 
disorder in periodically driven systems can induce AL in 
the time evolution. 

Controlling interactions between bosonic particles in 
the presence of a space periodic potential allows for in¬ 
vestigation of quantum phase transitions Q . In the limit 
of weak interactions, particles in the ground state reveal 
long-range phase coherence. For strong repulsive contact 
interactions, the phase coherence is lost because it is en¬ 
ergetically favorable to suppress quantum fluctuations of 
number of particles in each site of an external periodic 
potential. Such a Mott insulator (MI) phase is charac¬ 
terized by a gap in the excitation spectrum 0,0 . Tran¬ 
sition between the superfluid and Mott insulator phases 
has been demonstrated 0 and is extensively investigated 
in ultra-cold atom laboratories 0. In the present letter 
it is shown that repulsive particle interactions of a period¬ 
ically driven system lead to formation of Mott insulator¬ 
like state where long-time phase coherence is lost. 


process is currently a subject of the debate in the litera¬ 
ture In our study we do not consider the problem 

of the time crystal formation but concentrate on an anal¬ 
ysis of systems where time-periodicity is already given by 
an external driving. 


Results 


Anderson localization in the time domain. We 

consider single-particle systems that periodically depend 
on time, i.e. Hamiltonians Hit -\-T) = Hit) where T 
is the time period. While the energy is not conserved, 
there are so-called quasi-energy states (in analogy 
to Bloch states in spatially periodic problems) that are 
time-periodic eigenstates, HpUn = SnUn, of the so-called 
Floquet Hamiltonian Hpit) = H{t) — ihdt. There are 
many different quasi-energy eigenstates of the systems 
but we will be interested in those that are represented 
by wave-packets localized on classical s-resonant orbits, 
i.e. on orbits whose periods equal sT where s is inte¬ 
ger. Such extraordinary states exist in different experi¬ 
mentally attainable systems like, e.g., electronic motion 
in a hydrogen atom in microwave field rotating 

molecules [25| or an atom bouncing on an oscillating mir¬ 
ror The latter system will serve as an illustration 

of the ideas presented in this letter. In Fig. [1] classical 
motion of this system is described. From the semiclassi- 
cal point of view, existence of wave-packets which move 
along classical orbits is related to localization of quantum 
states inside elliptical resonant islands in the phase space 


It has been proposed that spontaneous breaking of 
time translation symmetry can lead to formation of time 
crystals where probability density at a fixed point in the 
configuration space reveals spontaneously time periodic 
behaviour |lll4l4l |. The possibility of such a spontaneous 


A single wave-packet moving on a s-resonant orbit can¬ 
not form a quasi-energy eigenstate because its period of 
motion is s times longer than the required period for all 
eigenstates of the Floquet Hamiltonian. However, su¬ 
perpositions of s wave-packets can form system eigen- 
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states - in Fig. [2] we show an example for the s = 4 case. 
There are s such superpositions that are linearly indepen¬ 
dent, hence, there are s eigenstates which reveal localized 
wave-packets moving along a s-resonant orbit. The cor¬ 
responding quasi-energy levels £r„ are nearly degenerated 
with only small splittings related to rates of tunneling of 
individual wave-packets. That is, if a single wave-packet 
is prepared initially on a s-resonant orbit, it travels along 
the orbit but after some time tunnels to the positions 
of other wave-packets which form the system eigenstates 
but which are missing initially. If s —>■ oo, the nearly de¬ 
generated eigenvalues e„ form a quasi-energy band which 
is an analogue of the lowest energy band of spatially peri¬ 
odic systems. Now we can establish conjecture about the 
behaviour of periodically driven systems and systems of 
particles in spatially periodic potentials. A single wave- 
packet localized on a s-resonant orbit is an analogue of 
a Wannier state localized in a single site of a spatially 
periodic potential [13 but in the time domain, see FigEl 
There exist also excited wave-packets that move on a s- 
resonant trajectory [53 which are analogues of Wannier 
states corresponding to excited energy bands of spatially 
periodic systems. The conjecture becomes formal when 
we derive Floquet energy of a periodically driven particle 
in the Hilbert subspace spanned by s individual wave- 
packets (j)j, 

Ef = ( 1 ) 

sT s 

where ((^|id_F|^)) = f and ■0 « ^ aj(j>j 

Q_ 

has been substituted [22]. The tunneling rates between 
wave-packets that are neighbours in the time domain (see 
Fig- mi Jj = ~‘^{{4>j+i\HF\4‘j)) with the same absolute 
value \Jj \ = J. 

Description of a resonantly driven particle has been 
reduced to the problem of a particle in ID lattice with 
nearest neighbour tunnelings ©. If we are able to create 

S 

an additional disorder term, E' = where Ej’s 

j=i 

are random numbers, Eq. ([1]) will become the ID Ander¬ 
son model and AL phenomena will emerge. The disorder 
term can be realized by an additional small perturbation 
E['(t) that fluctuates in time. That is, E['{t+sT) = E['{t) 
but between t and t + sT it behaves so that the set of 
Ej = {{(j)j\H'\cj)j)) reproduces a chosen set of random 
numbers. Then, the quasi-energy eigenstates of the total 
system, E[F{t) + are time periodic with the period 

sT but reveal superpositions of individual wave-packets 
(j)j with exponentially localized distributions. In Fig. [3] 
we describe an example of the realization of time disorder 
and show solutions that exhibit AL. 

Time crystal phenomena are related to time-periodic 
evolution of probability density for measurement of a par¬ 
ticle at a fixed position. If there is no time disorder, such 
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FIG. 1: Panel (a): schematic plot of a particle bouncing on 
an oscillating mirror in the presence of the gravitational field. 
In the absence of the oscillations all classical trajectories are 
periodic with a period the longer, the greater energy of a 
particle. When the time periodic oscillations are on and their 
amplitude is not very big, orbits whose period is multiple 
of the mirror oscillation period survive. Such perturbation 
resistant orbits are located in resonant islands in the classical 
phase space. If the period of an orbit is s times longer than the 
mirror oscillation period T, where s is integer, we will call the 
orbit s-resonant. Panel (b): description of the system is more 
convenient if we choose the coordinate frame that oscillates 
with the mirror [^ . Then, the mirror does not move but 
the gravitational acceleration becomes time dependent. The 
resulting Hamiltonian of the system reads H{t) = p^/2 -I- 
« -I- \zcos{2Ft/T) where the gravitational units have been 
used, i.e. Iq = (hd jvn?, to ~ and Eq = mglo 

for length, time and energy, respectively, where m is particle 
mass and g original gravitational acceleration. The results 
presented in Figs. 1 2131 correspond to A = 0.01 and T = 7.79. 


a probability reveals a uniform train of s humps that is 
repeated every sT period, see Fig. HKb). If the disor¬ 
der is on, the train is no longer uniform but reveals an 
Anderson localized distribution of humps, see Fig. |3l 
Mott insulator phase in the time domain. Let us 
consider N bosonic particles with repulsive contact in¬ 
teractions, characterized by a parameter go > 0, that 
are periodically driven, e.g. we can focus on ultra-cold 
atoms bouncing on an oscillating mirror. In the absence 
of the interactions, a Bose-Einstein condensate can be 
formed. Then, the system description reduces to a sin¬ 
gle particle problem and s-resonant driving can be de¬ 
scribed by the Floquet energy ([T]). In order to describe 
behaviour of the interacting many-body system we may 
truncate the Hilbert space to a subspace spanned by Fock 
states |ni,..., n^) where the occupied modes correspond 
to localized wave-packets 4>j moving along a s-resonant 
trajectory. Then, the many-body Floquet Hamiltonian 
reads 

1 ® 1 

ELf ~ 'y + h.c.) -f- — 'y ^ Uijfiifij, (2) 

where di and are bosonic anihilation and creation 
operators and fii = dldi. The coefficients Uij = 
goUYilYiYjlYj)) describe interactions between particles 
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FIG. 2: Analysis of the quasi-energy eigenstates localized on 
s-resonant periodic orbit, with s = 4, for a particle bouncing 
on an oscillating mirror - the classical behaviour is described 
in Fig. [T] Panel (a) shows one of the eigenstates at f = 0.25T 
(solid line) and t — O.ST (dash line). The eigenstate is a 
superposition of 4 wave-packets that move on the classical or¬ 
bit. Each wave-packet evolves with the period 4r but after 
T they exchange their positions what makes the whole eigen¬ 
state periodic with the period T. There are 4 such eigen¬ 
states that are localized on the orbit. Proper superpositions 
of the eigenstates allows one to extract 4 individual wave- 
packets, (j)j, that are numbered in (a) and (b). After each 
period T of the time evolution the wave-packets exchange 
their positions in the following order (fij+i (j>j (with j -I- 1 
modulo 4). Crystalline structures are not visible in the con¬ 
figuration space, however, they emerge in the time domain. 
Panel (b) shows time evolution of the 4 wave-packets (plotted 
with different colours), whose superpositions form the eigen¬ 
states localized on the 4-resonant orbit, at « = 121 that is 
close to the turning point of a particle in the classical descrip¬ 
tion. These wave-packets are analogues of Wannier states 
corresponding to the lowest energy band of a particle in a 
spatially periodic potential. Tunneling between wave-packets 
that are neighbours in the time domain are the leading tun¬ 
neling processes. These processes are taken into account in 
(0. The corresponding absolute values of the tunneling rates 
\Jj\ = 2\{{(j)j+i\HF\4>j))\ = J « 6 X 10“^ for the system pa¬ 
rameters described in Fig. 0 


that ocupy the same mode (for i = j) and between parti¬ 
cles in different modes {i ^ j). The latter are negligible 
if propagating wave-packets never overlap. If they pass 
each other during the evolution along the orbit, Uij ^ 0 
and the Hamiltonian (0) corresponds to a lattice model 
with long range interactions what could be surprising be¬ 
cause the original particle interactions are zero range. 
For the parameters described in Fig. 0 and s = 4, Uij are 
about O.lUii for i ^ j. 

If 0, the ground state of (0) is a superfluid state 
with long-time phase coherence. However, the long-time 
phase coherence is lost when Ua ^ NJ/s because the 



FIG. 3: Quantum description of a particle moving on the 
4-resonant orbit in the presence of a small perturbation, 
H'{t) = z otn cos {27r [n^ — sin (27r^)] }, that fluctu¬ 

ates in time. The whole Floquet Hamiltonian Hpit) -f H'(t) 
is time periodic with the period 4T and so do the eigen¬ 
states. The coefficients in H', are chosen so that the set 
of {{4>j\H'\4>3)) reproduces a chosen set of numbers Ej, where 
0j’s are the wave-packets described in Fig. 0 We have chosen 
Ej = 4Jwj where Wj are random numbers corresponding to 
a Lorentzian distribution W(wj) = l/(7r -|- TTWj). Then, the 
Floquet energy (0 supplemented with E' = Ej\aj\'^ con¬ 
stitutes the Lloyd model of ID lattice where all eigenstates 
are Anderson localized and the exact expression for the local¬ 
ization length is known Q. From the Lloyd model we know 
that the eigenstates are superpositions of the wave-packets, 
with |aj|^ oc where I is the AL length 

and jo is a number of the wave-packet around which a given 
eigenstate is localized. The wave-packets (j>j arrive at a given 
position 2 in equidistant intervals in time, thus, the AL length 
in time is It = IT. Solid lines in (a) and (b) show one of the 4 
eigenstates localized on the 4-resonant orbit at 2 = 121 versus 
time in the linear (a) and logarithmic (b) scales (the eigen¬ 
state is obtained in full numerical diagonallzation of Hf + H' 
but the solution of the Lloyd model is identical). Dash lines 
present behaviour of the eigenstates in the absence of the dis¬ 
order in time, i.e. when H' — 0. Despite the fact that the 
system is rather small, the characteristic exponential decay 
of the humps is clearly visible in (b) - the fitted exponential 
profile (dash-dotted line) corresponds to It = 0.18T. 


ground state becomes a Fock state |iV/s, N/s ,..., N/s) 
SE3- If we fix position in the configuration space and 
analyze what is the phase relation between wave-packets 
arriving one by one at this position (see Fig.0), it turns 
out the phase is totally undefined. Moreover, there is 
a gap in the excitation spectrum of the order of Ua. 
Hence, the system reveals Mott insulator-like properties 
but in the time domain. Ultra-cold dilute atomic gases 
are promising laboratories for the realization of such a 
MI phase due to an unprecedented level of experimental 
control. 
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FIG. 4: Comparison of crystalline structures in the space and 
time domains. Panel (a) illustrates a space periodic potential 
in ID with the periodic boundary conditions, the larger the 
system, the greater number of potential sites along the ring. 
Particles in the potential wells illustrate a Mott insulator state 
where in each potential site there is a well defined number 
of particles. Panel (b) refers to a crystal structure in the 
time domain with the time periodic conditions corresponding 
to the period sT. In the Mott insulator regime, at a fixed 
position, a well defined numbers of particles arrive every time 
interval T like on a conveyor belt or in a machine gun. The 
greater s, the larger time crystal is. 


Discussion 

Summarizing, we have shown that periodically driven 
systems can reveal non-trivial crystalline properties in 
the time domain. One-particle systems show Anderson 
localization in the time evolution if there are fluctuations 
in the periodic driving. Long-time phase coherence of pe¬ 
riodically driven many-body systems can be lost due to 
particle interactions and the Mott insulator-like phase 
emerges. 
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